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^The  III-V  semiconducting  alloys  are  typically  grown  by  epi¬ 
taxial  techniques  at  temperatures  where,  in  the  absence  of  sub¬ 
strate  effects,  they  are  thermodynamically  unstable.  This  can 
result  in  problems  associated  with  clustering  of  like  atoms  or 
ordering  of  unlike  atoms.  Long-range  ordering  could  yield 
interesting  III-V  ternary  compounds. 

The  mixing  enthalpy  of  III-V  semiconductor  alloys  is  fairly 
well  described  by  regular  solution  theory,  with  a  thermodynamic 
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interaction  parameter  that  is  sensitive  to  the  lattice  spacing  of 
the  binary  constituents.  An  estimate  of  the  interaction  parameter 
is  derived  from  a  model  which  ascribes  the  mixing  enthalpy  to  bond 
distortions  associated  with  the  alloy  formation,  and  relates  these 
to  the  macroscopic  elastic  properties  of  the  crystal.  Numerical 
estimates  are  given  for  the  18  alloys  with  cations  Al,  Ga,  In  and 
anions  P,  As,  Sb,  and  these  are  compared  with  experimental  values 
and  alternative  models.  To  within  a  single  adjustable  parameter. 


A  calculation  of  the  elastic  energy  associated  with  composition 
fluctuations  (clustering)  in  these  alloys  indicates  that  this  energy 
is  sufficient  to  suppress  clustering  above  the  critical  mixing  tem¬ 
perature.  For  semiconductor  alloys  AxBi_xC  where  the  lattice  para¬ 
meters  of  AC  and  BC  differ,  clustering  is  reduced  below  that  exhibited 
by  a  purely  random  distribution.  This  provides  a  tendency  toward 
order  in  these  alloys.  The  extent  to  which  clustering  is  reduced 
is  estimated  and  related  to  recent  observations  of  Raman  linewidths 
and  lineshapes  in  Ini_xGaxAs  and  Ali_xGaxAs. 

Below  bandgap  photovoltaic  response  with  a  half-power  point  at 
1.31eV  is  observed  in  p-type  Ini_xGaxP  on  n-type  GaAs  he tero junctions 
with  0.11}  interfaces.  This  response  is  apparently  due  to  photon- 
assisted  tunneling  of  carriers  across  the  interfacial  energy  gap, 
which  is  smaller  than  the  energy  gaps  of  the  constituent  materials. 

The  same  mechanism  can  be  employed  in  other  he tero junctions  with 
appropriate  energy  band  lineups  to  obtain  longer  wavelength  infrared 
response. 
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1.  INTRODUCTION 


Although  there  is  substantial  interest  in  III-V  semi¬ 
conductor  alloys  for  electronic  and  optoelectronic  devices, 
relatively  little  is  known  about  the  distribution  of  atoms 
in  these  materials  or  the  effects  of  non-random  distributions 
on  device  performance.  In  many  of  these  alloys,  however, 
non-random  distributions  of  atoms  are  expected  to  be  important 
from  thermodynamic  considerations.  This  is  particularly  true 
at  the  low  temperatures  commonly  used  for  epitaxial  growth. 

The  available  evidence  suggests  that  short-range  clustering 
of  like  atoms  or  short-range  ordering  of  unlike  atoms  can 
produce  device  problems  such  as  excess  noise  and  leakage 
current,  premature  voltage  breakdown,  and  lower  carrier 
mobilities.  Long-range  ordering  of  unlike  atoms,  however, 
could  potentially  yield  III-V  ternary  compounds  with  proper¬ 
ties  superior  to  their  parent  alloys.  In  epitaxial  growth, 
substrate  effects  such  as  lattice  match  and  nonequivalent 
sublattice  sites  are  expected  to  have  a  strong  influence  on 
these  phenomena.  Also,  two  of  the  important  III-V  alloys, 
GaJ{In1_xP  and  GaxIni_xAs/  can  be  epitaxially  grown  lattice 
matched  to  GaAs  and  InP,  respectively,  at  compositions  near 
the  point  of  greatest  ordering  probability.  The  objective 
of  this  work  is  to  investigate  various  aspects  of  clustering 
and  ordering  in  these  materials. 


2.  ELASTIC  MODEL  OF  III-V  ALLOY  FORMATION 

The  calculation  of  phase  diagrams  and  of  other  thermodynamic 
properties  of  tetrahedral  ternary  III-V  semiconductor  alloys 
of  the  form  AxB^_xC  (or  ACxD^_x>  using  regular  solution  theory 
has  proved  to  be  a  useful  and  reasonably  accurate  approximation. 
Several  of  the  ternary  alloys  with  cations  Al,  Ga,  In  and  anions 
P,  As,  Sb  have  important  present  or  potential  applications  in 
electronic  technology.  These  eighteen  alloys  are  mixtures  of 
binary  compounds  with,  in  general,  different  lattice  parameters. 
The  lattice  parameters  of  the  binary  compounds  are  different 
when  the  AC  and  BC  bond  lengths  are  different,  and  the  lattice 
parameter  of  the  alloy  then  differs  from  both.  The  lattice  para¬ 
meters  of  all  the  ternary  alloys  that  can  be  formed  from  the  ions 
listed  except  those  involving  the  cation  pair  Al,  Ga  are  sensi¬ 
tive  to  composition.  One  then  may  take  the  view  that  the  bonds 
are  distorted  or  rearranged  in  forming  the  alloy.  The  energy 
associated  with  the  bond  distortion  can  be  estimated  from  the 
macroscopic  elastic  properties  of  the  crystals.  We  discuss 
the  contribution  of  this  energy  to  the  enthalpy  of  mixing  and 
its  effect  on  clustering  and  spinodal  decomposition.  We  suggest 
that  it  is  a  major  or  dominant  factor  in  both  phenomena. 

2 . 1  MIXING  ENTHALPY 

In  regular  solution  theory  the  molar  enthalpy  of  mixing  of 


a  binary  solution  is  given  by 


-3- 


where  fl  is  the  interaction  parameter  and  where  x  and  1  -  x  are 
the  fractions  of  the  components  making  up  the  solution.  In 
terms  of  a  picture  of  pairwise  interaction  energies,  the  inter¬ 
action  parameter  is 

ft  =  6NQe  (2) 

where  6Nq  is  the  number  of  interacting  pairs  on  the  fee  sublat¬ 
tice  and 

6  ■  2eAB  '  £AA  -  eBB  (3) 

is  the  excess  of  the  AA  and  BB  bond  energies  between  like  ions 
over  the  AB  bond  energy  of  unlike  ions.  In  the  ternary  alloys 
under  discussion  these  second  nearest  neighbor  interactions  are 
thought  to  be  rather  small,  consistent  with  the  small  experi¬ 
mental  values  of  the  interaction  parameter. 

Stringfellow,  using  the  covalent  bonding  ideas  of  Phillips 

and  Van  Vechten1^ ,  has  introduced  the  delta  lattice  parameter 
2) 

model  which  provides  values  of  the  interaction  parameter  of 
all  the  alloys  in  terms  of  a  single  adjustable  parameter.  The 
agreement  between  the  DLP  theory  and  experiment  is  remarkably 
good  for  a  theory  with  a  single  adjustable  parameter.  The 
model  takes  a  less  specific  approach  than  the  identification 
of  Equations  (2)  and  (3)  of  the  interaction  parameter  with  the 
sum  of  binary  bond  energies,  and  assumes  that  the  energy  is  a 
more  global  property  of  the  semiconductor  that  depends  only  on 
the  lattice  parameter  of  the  crystal.  The  idea  is  that  the 
bonding  energy  is  proportional  to  the  band  gap  which,  in  turn. 


is  proportional  to  a  power  of  the  lattice  parameter.  This  is 

3\ 

similar  to  but  distinct  from  Van  Vechten's  proposal  '  that  is 
produced  by  the  reduction  of  the  band  gaps  due  to  the  disorder 
associated  with  alloying.  The  DLP  model  leads  to  an  approxi¬ 
mate  value  for  fi  of 

fi  =  <35/8)KA2d  "4‘5  (4) 

s 

where  A  =  a  -  b,  5  =  1/2 (a+b) ,  with  a  and  b  the  lattice 
parameters  of  AC  and  BC,  respectively.  The  quantity  K  is  ob¬ 
tained  by  a  least  squares  fit  to  the  experimental  data  and  is 

7  °2  5 

equal  to  1.15  x io  cal/mole  A  *  .  The  excellent  fit  to  the 

2 

experimental  data  depends  largely  on  the  A  dependence  of  ft. 

We  wish  to  show  here  that  if  part  or  all  of  ft  is  inter¬ 
preted  as  being  associated  with  strain  or  lattice  parameter 

mismatch,  then  it  is  readily  shown  that  this  contribution 

2 

exhibits  the  same  A  dependence. 

It  is  observed  that  the  average  lattice  spacing  of  the 
alloys  closely  follows  Vegard's  law 

dQ  =  xa  +  ( 1  -  x)  b  .  (5) 

Since  the  difference  between  the  lattice  parameters  a  and  b  of 
the  binary  compounds  corresponds  to  the  different  AC  and  BC 
bond  lengths,  the  bonds  must  be  distorted  or  rearranged  to  form 
the  alloy.  If  a  and  b  do  not  differ  too  widely,  the  energy  of 
the  rearranged  configuration  must  be  of  second  order  in  a  -  dQ 
and  b  -  d  .  Thus  the  contribution  to  the  enthalpy  of  formation 


of  the  alloy  is 


AH  =  (x(d  -  a)  2  +  (1  -  x)  (d  -  b)  2)  (6) 

in  ^  o  o 

where  X  is  an  adjustable  parameter.  With  dQ  from  Equation  (5) 
this  becomes 


AH  -  x(l-x)A  (|)2  (7) 

m  a 

2 

with  the  same  A  dependence  as  the  DLP  model. 

Since  we  are  viewing  the  effect  as  an  elastic  distortion, 
we  can  obtain  a  crude  numerical  estimate  for  X  from  tne  macro¬ 


scopic  elastic  properties  of  the  crystal.  At  one  extreme  one 
can  assume  a  strict  virtual  crystal  model  of  the  alloy,  which 
would  identify  the  bond  distortion  energy  with  that  of  an 
isotropic  compression  or  dilatation.  At  x  *  1/2,  the  principal 
strains  are  given  by  1/2A  and  the  elastic  energy  per  unit 


volume  is 


yielding 


W  =  1/2B  (1  A)  2 
z  d 


(8) 


X  «  i  BV,  Q*  =  A  (4)  2  (9) 

*  d 

where  B  is  the  bulk  modulus  and  V  is  the  molar  volume  of  the 

alloy  with  x  =  1/2.  Table  2.1  lists  ft  and  ft '  for  the  ternary 

s 

compounds  under  discussion. 

Equation  (9),  which  contains  no  adjustable  parameters, 
yields  higher  interaction  parameter  values  than  the  DLP  theory. 
A  much  better  fit  to  the  experimental  data  is  obtained  if  the 


prediction  of  Equation  (9)  is  reduced  by  a  factor  of  about 

4.4,  resulting  in  the  column  labeled  ft' '  in  Table  2.1.  Indeed, 

the  entries  in  that  column  fit  the  observed  values  as  closely 

2) 

as  Stringfellow' s  "best  fit"  2.45  power  law.  Moreover,  a 
plausible  case  can  be  argued  for  a  reduction  of  the  estimate 
by  a  factor  of  this  order. 

4) 

It  has  been  demonstrated  recently  that  although  the 
average  lattice  parameter  closely  follows  Vegard's  law,  the 
AC  and  BC  bond  lengths  do  not  take  the  common  nearest  neighbor 
value  of  the  virtual  crystal  but  remain  much  closer,  by  a  factor 
of  4  or  5,  to  their  lengths  in  the  binary  compounds.  If  it  is 
assumed  that  the  stretching  and  compression  of  nearest  neighbor 
bonds  is  the  only  distortion,  this  reduces  the  estimate  of  A.  by 
the  square  of  this  factor,  of  order  20.  This,  however,  almost 
certainly  underestimates  the  effect  of  the  distortion,  since  the 
conservation  of  nearest  neighbor  bond  lengths  requires  a  dis¬ 
tortion  of  the  bond  angles  and  of  second  neighbor  bond  lengths 
and  angles. 

In  addition,  the  fact  that  the  AC  and  BC  bond  lengths  tend 
to  remain  constant  (to  about  80%)  independent  of  alloy  compo¬ 
sition,  necessarily  implies  that  the  A  and  B  atoms  in  the  alloy 
are  much  more  correlated  than  in  a  random  solution.  These  strong 
correlations  could  well  be  due  to  the  strong  elastic  forces. 
Neither  our  theory  nor  any  other  theory  currently  takes  this 
into  account  properly.  Thus  within  the  scope  of  our  model, 


Table  2.1  Experimental  and  calculated  values  of  the  interaction 

parameter.  Bulk  moduli  by  interpolation  between  binary 
compound  values  from  References  8  and  9;  interaction 
parameters  (in  k  ca 1/mole)  in  Columns  3,  4,  and  5  quo¬ 
ted  from  Stringf ellow,  Reference  2  (small  differences 
in  the  calculated  values  in  Column  4  because  we  use 
fewer  significant  numbers  than  Reference  2  for  A/d) . 


Table  2.1  furnishes  strong  evidence  that  this  quasi-elastic 
energy  is  a  significant,  and  perhaps  a  dominant  part  of  the 
enthalpy  of  mixing. 

2.2  COMPOSITION  FLUCTUATIONS 

A  binary  mixture  whose  enthalpy  of  mixing  is  positive  has 
a  tendency  to  separate  into  its  components.  This  tendency  mani¬ 
fests  itself  in  the  form  of  spinodal  decomposition  and  the 
development  of  a  miscibility  gap  below  a  critical  mixing  tem¬ 
perature  and  in  the  form  of  clustering  on  a  microscopic  scale 
above  this  temperature.  In  the  alloys  under  discussion  here, 
this  tendency  is  most  readily  interpreted  in  terms  of  prefer¬ 
ential  pairwise  interactions  as  expressed  by  Equations  (2)  and 
(3) ,  but  it  does  not  depend  on  this  interpretation.  Jones, 
Porod,  and  Ferry5^  have  recently  examined  the  degree  of  cluster¬ 
ing  that  may  be  expected  to  occur  above  the  critical  temperature 
of  mixing  as  a  function  of  the  interaction  parameter  interpreted 
in  this  manner,  of  the  temperature,  and  of  the  mole  fraction  x. 
We  extend  this  calculation  to  include  an  elastic  energy  which  is 
large  enough  to  overcome  the  clustering  tendency. 

Since  the  lattice  parameters  of  the  alloys  follow  Vegard's 
law,  the  ions  of  the  alloy  at  least  on  the  average  occupy  po¬ 
sitions  on  a  virtual  crystal  lattice  with  the  lattice  parameter 

of  Equation  (5).  Note  that  this  statement  is  not  in  conflict 

4) 

with  the  observation  cited  above  that  the  bond  lengths  need 
not  follow  Vegard's  law. 


Clustering  produces  a  fluctuation  from  the  average  compo¬ 
sition,  while  the  region  of  the  cluster  must  occupy  the  volume 
prescribed  by  the  average  virtual  crystal  lattice.  Thus  the 
region  of  the  cluster  is  strained.  The  energy  associated  with 
the  strain  tends  to  reduce  the  degree  of  clustering,  but  we 
shall  see  that  it  does  not  eliminate  a  miscibility  gap  below 

the  critical  mixing  temperature.  The  effect  is  similar  to  the 

6 )  7) 

"lattice  latching"  and  stabilization  against  decomposition 

associated  with  the  mismatch  strains  of  heteroepitaxy.  In  the 
present  instance  we  deal  with  the  strains  arising  from  fluctua¬ 
tions  in  a  homogeneous  medium.  In  either  case,  the  existence 
of  a  strained  region  presupposes  the  existence  of  a  crystal 
lattice  to  be  matched  by  the  continuing  growth,  in  epitaxy, 
this  lattice  is  furnished  by  the  substrate.  In  homogeneous 
growth  it  depends  on  the  formation  of  a  nucleus,  which  at  a 
given  composition  can  only  form  at  temperatures  above  the 
spinodal  decomposition  curve.  This  is  the  reason  why  the  strain 
energy  affects  clustering  but  not  the  miscibility  gap. 

Let  the  average  composition  of  an  alloy  be  A  B.  C  and 

xo  -L“x0 

the  local  composition  AxBi_xc-  Then  the  lattice  parameter  is 

d  =  x  a  +  (1  -  x  ) b.  (10) 

o  o  o 

If  the  crystal  were  permitted  to  relax  everywhere,  the  local 
lattice  parameter  would  be 


However/  the  compatibility  constraint  imposes  a  local  strain 

6  =  d  -  dQ  =  (x  -  xQ)  {a  -  b)  =  £  A  ,  (12) 

the  second  equality  defining  the  composition  fluctuation  £. 

The  energy  density  of  the  strain  can  be  computed  as  in 
Equation  (8)  and  is 

W  =  (9/2)  B£2  2  (13) 

d 

per  unit  volume. 

Unlike  Equation  (8)  for  the  elastic  contribution  to  the 
mixing  enthalpy.  Equation  (12)  should  not  be  a  significant 
overestimate  of  the  strain  energy.  In  the  present  application 
we  are  computing  the  energy  of  actual  compressions  and  dilata¬ 
tions  in  terms  of  measured  elastic  constants  and  the  nature  of 
the  bond  distortions  accompanying  the  strain  is  irrelevant. 
Therefore,  if  the  mixing  enthalpy  is  dominated  by  the  bond 
distortion  effects  discussed  in  the  last  section,  then  it  is 
evident  that  the  strain  energy  of  composition  fluctuations  is 
more  than  adequate  to  overcome  the  clustering  tendency. 

A  similar  argument  can  be  made  if  the  clustering  is  ascribed 
to  preferential  pairwise  interactions.  The  molar  bond  energy 
density  of  the  ternary  alloy,  taking  into  account  the  4  nearest 
and  12  next  nearest  neighbor  interactions  is 

E  =  NQ(4eAx  +  4eB (1  -  x)  +  Se^x2  +  12eABx(l  -  x)+6efiB (1  -  x)  2) ,  (1 

where  eA  and  e0  are  the  AC  and  BC  interaction  energies. 


The  contribution  to  the  bond  energy  density  of  a  composition 
fluctuation  £  is  E(xQ+£)  -E(xQ)  =  n,  given  by 

n  =  N0((4eA  -  4efi  +  126^  -  12eBB  -  12  xQ  e)  £  -  6e£2)  (15) 

with  e  as  defined  in  Equation  (3) . 

In  a  crystal  of  a  given  composition,  the  contribution  to 

the  interaction  energy  linear  in  the  fluctuations  must  by  defi- 

2 

nition  average  to  zero.  Therefore  only  the  term  in  £  need  be 
considered  for  the  purpose  of  minimizing  the  free  energy.  But 
from  Equation  (2)  this  term  is  just 

n  -  -  nz,2.  (16) 

so  that  the  net  molar  energy  of  the  fluctuations,  from 
Equations  (13)  and  (16),  is 

wv  +  n  -  (o'  -  n) £2  (17) 

averaged  over  the  crystal. 

Equation  (17)  yields  a  positive  (or  zero)  fluctuation 
energy  for  each  of  the  ternary  alloys  we  are  considering,  indi¬ 
cating  that  clustering  is  suppressed  by  the  strain  energy. 

2.3  CONCLUSIONS 

In  our  model  the  formation  of  a  tetrahedral  semiconductor 
alloy  is  accompanied  by  bond  distortion,  and  an  estimate  of  the 
electronic  energy  of  this  distortion  can  be  obtained  from  the 
macroscopic  elastic  properties  of  the  crystal.  Superficially 
this  view  appears  distinct  from  the  basis  of  the  DLP  model, 
which  is  deduced  from  the  empirical  lattice  parameter  dependence 


f 


of  the  optical  bandgap  and  from  Van  Vech ten's  proposal.  But 
the  bandgap  itself,  as  pointed  out  by  Phillips  and  Van  Vechten^ 
is  a  measure  of  the  covalent  bonding  energy.  The  elastic  constants 
of  a  crystal  depend  as  directly  on  the  bonding  energy  as  its  opti¬ 
cal  and  dielectric  properties.  Therefore  it  is  perhaps  not  too 
surprising  that  models  derived  from  optical  and  from  elastic  crys¬ 
tal  parameters  should  yield  predictions  of  similar  form  for  the 
thermodynamic  properties.  The  merit  we  may  claim  for  the  elastic 
model  is  that  it  contributes  a  new  and  convenient  experimental 
approach  to  the  thermodynamic  quantities  and  that  it  may  yield 
useful  theoretical  insights. 
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3.  COMPOSITION  CORRELATIONS  IN  ALLOYS 


Ternary  tetrahedrally  coordinated  semiconductor  alloys 
AxB^_xC  or  ACxDi_x  where  A  and  B  are  group  III  (or  II)  ele¬ 
ments  and  C  and  D  group  V  (or  VI)  elements  can  for  many  pur¬ 
poses  be  described  as  binary  mixtures  whose  component  particles 
reside  (approximately)  on  the  sites  of  an  fee  lattice. 

Specifically,  thermodynamic  parameters  such  as  mixing  enthalpy ^ 

(2) 

and  spinodal  curves  conform  fairly  well  to  the  rules  de¬ 
scribing  regular  solutions. 

The  key  quantity  of  regular  solution  theory  is  the  inter¬ 
action  parameter  ft,  which  can  be  interpreted  in  terms  of  pair 
interaction  energies,  thus 

V(2eab  -  eaa  -  EBB>' 

where  Nq  is  Avogadro's  number,  Z  is  the  coordination  number 

(12  for  fee)  and  E. .  are  the  (negative)  pair  interaction 

1 J 

energies.  A  negative  value  of  ft  indicates  a  tendency  to 
ordering  (compound  formation) ,  a  positive  value  a  tendency 
to  decomposition.  All  the  observed  interaction  parameters 
of  the  alloy  group  with  A  =  A1,  Ga,  In;  B  *  P,  As,  Sb; 

A=Zn,  Cd,  Hg;  B  =  Se,  Te  are  positive  or  zero.  Thus  one 
expects  -  and  observes  -  miscibility  gaps  and  spinodal  de¬ 
composition  below  a  critical  mixing  temperature  characteristic 
of  each  alloy. 

The  observed  interaction  parameters  are  found  to  be 
proportional  to  the  square  of  the  lattice  parameter  difference 
between  the  pure  compounds  that  constitute  the  alloy.  This 


dependence  has  been  interpreted,  in  the  so-called  delta- 
lattice-parameter  (DLP)  model  in  terms  of  the  lattice  para¬ 
meter  dependence  of  the  average  bandgap^,  and  more  recently, 
especially  after  details  of  the  composition  dependence  of 

bond  lengths  in  Ga^_xInxAs  were  revealed  by  EXAFS  measure- 

(4) 

ments  ,  in  terms  of  the  bond  distortions  associated  with 
the  formation  of  the  alloy.  By  relating  the  bond  distortion 
energies  to  the  macroscopic  elastic  constants  of  the  crystals, 
a  good  case  could  be  made  for  ascribing  most  or  all  of  the 
observed  mixing  enthalpy  to  this  cause^'^. 

A  positive  mixing  enthalpy  would  normally  produce  cluster 

ing  of  like  particles  above  the  critical  mixing  temperature, 

and  an  estimate  of  the  degree  of  clustering  to  be  expected 

(7) 

was  published  recently  .  For  the  model  that  was  used  in 
this  calculation,  the  origin  of  the  positive  interaction  para¬ 
meter  is  irrelevant.  It  has  been  pointed  out,  however,  that 
in  the  type  of  "lattice  gas"  that  represents  these  alloys, 
with  a  composition  dependent  lattice  parameter,  if  an  alloy 
crystal  is  formed  at  all,  elastic  interactions  will  suppress 
the  clustering ^ .  It  was  also  suggested  that  these  elastic 
interactions  were  likely  to  produce  an  ordering  tendency, 
that  is  to  say,  that  local  composition  fluctuations  might  be 
reduced  below  the  level  to  be  expected  of  a  purely  random 
distribution  of  the  constituents.  What  follows  is  a  simple 
thermodynamic  model  for  the  composition  correlations  induced 
by  the  elastic  energy. 
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3 . 1  THERMODYNAMIC  MODEL 

The  lattice  parameters  of  these  alloys  follow  Vegard's 

law  quite  accurately.  Any  deviation  from  a  linear  dependence 
of  lattice  parameter  on  composition  (bowing)  is  small  enough 
to  be  quite  negligible  for  the  purpose  of  this  argument. 

Consider  a  region  of  the  alloy  large  enough  to  be  regarded 
as  macroscopic.  If  the  composition  of  this  region  deviates 
from  the  average  composition  of  the  alloy,  but  must  be  accomo¬ 
dated  in  the  volume  of  the  virtual  lattice  that  "belongs  to 
it",  it  will  suffer  compression  or  dilatation.  This  is  illustra^ 
ted  schematically  in  Figure  3.1. 

The  energy  associated  with  this  strain  is  readily  computed*5 

Let  the  average  composition  of  the  alloy  be  A  B.  C  and  the 

o  o 

local  composition  A  B.  C.  Then  the  lattice  parameter  is 

dQ  =  xQ  a  +  (l-xQ)b  (1) 

where  a  is  the  lattice  parameter  of  AC  and  b  the  lattice  para¬ 
meter  of  BC.  If  the  crystal  were  permitted  to  relax  every¬ 
where,  the  local  lattice  parameter  would  be 

d  *  xa  +  (l-x)b,  (2) 

however,  the  compatibility  constraint  imposes  a  local 
strain 


6  =  d-d^  = 
o 


(x-x_J (a-b)  3  Ax  (a-b) 


(3) 
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wit  h  an  associated  energy  density  (per  unit  volume) 


W( Ax) 2  =  18B  (Jj|)  (Ax)2 


(4) 


where  B  is  the  bulk  modulus  of  the  crystal. 

Equation  (4)  is  based  on  isotropic  compression  or 

dilation;  an  analogous  expression  can  be  obtained  for  a 

planar  strain.  In  that  form  the  elastic  energy  is  thought 

to  be  responsible  for  the  "lattice  latching"  observed  in 
(8) 

heteroepitaxy  ,  the  growth  of  a  single  lattice- 
matched  alloy  composition  over  a  range  of  melt  compositions. 
It  also  accounts  for  the  successful  heteroepitaxial  growth 
of  alloy  compositions  unstable  in  the  bulk  at  the  growth 
temperature . 

It  has  been  showiPlhat  a  reasonable  fit  to  the  inter¬ 
action  parameters  of  the  III-V  alloys  considered  here  is 
given  by  the  semi-empirical  formula 


W 


0.226  X  18B  ( 


a-bN 

a+K 


(5) 


making  a  contribution  the  the  mixing  enthalpy  per  mole  of 


AH  -  W  x  ( 1— x)  V 

where  V  is  the  molar  volume  of  the  alloy. 

The  entropy  of  mixing  is  approximately 


S  =  In 


V  NB! 


naa1  nbb!  (nab!) 


(6) 


(7) 


where  the  are  numbers  of  ions  and  the  numbers  of 
second  neighbor  pairs.  Per  mole  this  is 

S  *  R[xlnx  +  (l-x)ln(l-x)  -x2ln(x2)  -  (1-x)  2ln(l-x)  2 

-  2x ( 1-x) lnx ( 1-x) ] .  ( 8 ) 

Using  Equations  (4),  (6),  and  (8)  in  the  thermodynamic  for- 
(9) 

mula  for  macroscopic  fluctuations  of  concentration  on  an 
assembly  of  N  sites  with  total  Gibbs  free  energy  G 

<  (Ax)  2>  =  -  (8) 

o2  c/ax2)T^N 

we  find 


<  (Ax) 


"  n  X  (1~X) 


U  + 


2 (W-W  )V 
RT 


x(l-x) ] 


-1 


(10a 


or,  in  terms  of  the  critical  mixing  temperature  Tc 
given  by  2RT  =w'v  =  0 . 226WV 

C 


<(Ax)2>=i 


x(l-x) 


-1 


|l  +  13.7^x(l-x2J 


(10b 


One  may  guess  conservatively  that  the  disorder  grown  into 
a  crystal  corresponds  roughly  to  a  value  frozen  in  at  the 
growth  temperature.  Thus  if  we  compare  a  lattice-matched 
alloy  such  as  GaQ  jAIq  ^As  with  Tc  =  0  with  the  mismatched 

O 

alloy  Ing  5As  with  Tc  ~  600  K  grown  epitaxially  at 

O 

T  -  900  K,  we  expect,  from  Equation  (10a)  or  (10b)  to  find 


a  relative  reduction  of  the  mean  square  concentration 


fluctuations  in  the  mismatched  alloy  by  a  factor  of  ~  0.3. 

The  reduction  would  be  expected  to  be  more  significant, 
depending  on  cooling  rate  and  on  other  growth  parameters, 
if  appreciable  solid  state  diffusion  takes  place  after  growth. 

3.2  RAMAN  SCATTERING  ANALYSIS 

Some  qualitative  experimental  support  for  this  prediction 
is  to  be  found  in  linewidths  and  lineshapes  observed  by  Parayanthal 
and  Poliak *10^  in  Raman  scattering  from  A1xGai_xAs  and  InJCGa1  As. 
In  this  work  it  was  found  that  the  broadening  and  asymmetry  of 
the  Raman  line  was  more  pronounced  in  the  lattice-matched  alloy 
A1  Ga,  As  than  in  the  mismatched  In  Ga.  As. 

X  1“X  X  X-X 

A  "spatial  correlation"  model  has  been  successfully  used 

to  interpret  the  effect  of  microcrystalline  ,  and  implanta- 
(12) 

tion  damage'  disorder  on  Raman  linewidths  and  lineshapes. 

This  model  replaces  the  phonons  of  the  infinite  ordered 
crystal  with  phonons  localized  within  "correlation  regions" 
whose  size  is  a  measure  of  the  partial  ordering  in  the  imperfect 
crystal.  As  a  result,  the  q  =  0  momentum  selection  rule  of 
the  perfect  crystal  is  relaxed,  and  phonons  in  a  region  of  the 
Brillouin  zone  corresponding  to  the  correlation  length  par¬ 
ticipate  in  the  scattering.  Experimental  values  of  the  cor¬ 
relation  length  are  deduced  by  fitting  the  observed  linewidths 
and  lineshapes,  assuming  a  gaussian  size  distribution  of  the 
correlation  regions,  generally  in  reasonable  agreement  with 
the  estimated  sizes  of  crystallites^11^  or  undamaged  regions^12^. 


Substitutional  disorder  in  an  alloy  similarly  breaks  the 
crystal's  translational  symmetry  and  hence  the  q  =  0  selection 
rule  for  Raman  scattering.  Since  the  perturbation  is  the 
phonon  modulation  of  the  dielectric  tensor,  the  effective  size 
of  the  scattering  regions  is  given  by  the  correlation  length 
of  this  tensor's  spatial  fluctuations  .  No  theoretical 
prediction  has  so  far  been  made  for  this  quantity,  but 
Parayanthal  and  Poliak  were  able  to  fit  their  observed 
Raman  lines  with  a  spatial  correlation  model  in  which  the 
correlation  length  is  simply  a  phenomenological  parameter 
determined  from  the  experimental  data.  For  their  epitaxially 
grown  alloys,  this  parameter  ranges  from  about  8  to  above  30 
lattice  parameters. 

Of  the  two  alloy  series  for  which  results  are  given, 
GaJCAl1_xAs  has,  within  experimental  accuracy,  zero  mixing 
enthalpy,  and  the  cations  should  constitute  a  perfect  lattice 
gas  on  their  fee  sublattice.  Their  concentration  fluctuations, 
which  may  be  taken  as  a  measure  of  the  disorder,  then  are 
given  by  Equations  (10)  with  Tc=0.  The  Raman  lines  bear 
this  out  qualitatively,  in  the  sense  that  the  linewidths 
and  asymmetries  are  greatest  for  samples  with  x  near  0.5. 

The  broadening  of  the  Raman  lines  even  in  the  completely 
random  alloy  Ga  Al.  As  is  modest  as  compared  with  the 
spectrum  of  amorphous  materials,  and  the  spatial  correlation 
model  is  successful  in  the  interpretation  of  the  lineshapes. 
This  suggests  that  one  might  assign  the  role  of  the  defects 


that  disrupt  the  phonon  modes  to  fairly  large  fluctuations 
from  the  average  composition.  The  correlation  or  coherence 
length  can  then  be  related  to  the  statistics  of  large 
fluctuations. 

In  a  random  binary  alloy,  the  probability  distribution 
of  ions  is  binomial.  With  a  coordination  number  Z  =  12,  the 
distribution  of  even  the  nearest  neighbor  shell  is  fairly 
well  represented  by  a  Gaussian,  and  the  probability  of 
occurrence  of  a  deviation  1  5  is  approximated  by 

P<£)  =  erfc  e  (11) 

where  £  =  £//  2  < (Ax)  >  . 

In  order  to  arrive  at  an  order  of  magnitude  estimate 
of  the  correlation  length,  we  take  the  size  of  the  defect  to 
comprise  one  cation  shell  of  Z  ions.  If  a  deviation  of 
relative  magnitude  £  is  to  occur  in  a  region  containing 
NZ  cations,  we  should  have  NZP(£)  =  1,  and  therefore  the 
average  distance  between  such  defects,  which  we  take  to  be 
the  correlation  length,  is 

L  =  N1/3a=  (ZerfcC)"1/3a  (12) 

where  a  is  the  lattice  parameter.  Equation  (12)  yields 
values  consistent  with  the  estimates  based  on  the  Raman 
spectra  if  £  is  about  two  standard  deviations.  Using  the 
average  value  of  L  for  the  GaAs-like  mode  in  Ga  Al^  As  for 
the  four  samples  with  x  near  0.5,  we  find  5  =  1.88. 


Equations  (10)  and  (12)  also  provide  a  comparison  of  the 


correlation  lengths  in  Ga  In.  As  and  GaAl.  As.  Because 

X  1“X  X  JL  —  X 

of  the  suppression  of  concentration  fluctuations,  we  expect 
the  scattering  centers  associated  with  wide  fluctuations  to 
be  rarer,  increasing  the  correlation  length.  Using  the  fac¬ 
tor  0.3  estimated  from  Equation  (10)  for  the  reduction  in  the 
mean  square  fluctuation  in  Equation  (12) ,  we  find 


Ga0.5In0.5As 

JGa0.5A10.5As 


«  12. 


(13) 


This  factor  is  larger  than  the  numbers  deduced  from  the  Raman 
spectra  by  Parayanthal  and  Poliak,  suggesting  that  additional 
phonon  scattering  processes  are  active. 

3.3  X-RAY  LAUE  PATTERNS 


The  reduction  in  randomness  discussed  above  is  expected  to 
have  other  experimentally  testable  consequences.  For  example, 
it  is  possible  to  examine  such  effects  with  x-ray  back-reflection 
Laue  photographs,  where  long-range  order  in  a  III-V  alloy  would 
produce  reflections  other  than  those  expected  for  the  sphalerite 
crystal  structure.  Figure  3.2(a)  shows  a  Laue  pattern  for  an 
AIq  4GaQ  gAs  epitaxial  layer  grown  on  a  GaAs  substrate,  while 
Figure  3.2(b)  shows  a  pattern  with  the  same  crystal  orientation 
for  an  InQ  5Ga0  epitaxial  layer.  No  additional  spots  are  ob¬ 
served  for  the  InQ  gP  over  the  AlQ  4GaQ  gAs,  indicating  no 

long-range  order.  The  064  and  033  reflections,  however,  are 
much  more  intense  for  InQ  5GaQ  gP  than  for  AlQ  5GaQ  5As*  We 
are  currently  examining  the  implications  of  these  measurements 
for  correlations  and  ordering  in  these  alloys. 
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BELOW  BANDGAP  PHOTORESPONSE  IN 


Ga„P-GaAs  HETEROJUNCTIONS 


Several  schemes  have  been  used  in  intrinsic  infrared 

photodetectors  to  obtain  response  below  the  bandgap  of  the 

constituent  semiconductors.  One  of  these  methods  is  to  form 

an  alloy  between  two  compounds  which  have  conduction  and 

valence  bands  of  opposite  symmetry,  such  as  PbTe  and  SnTe*. 

Since  the  conduction  and  valence  bands  must  interchange  at 

some  alloy  composition,  they  move  toward  one  another  producing 

alloy  bandgaps  smaller  than  the  endpoints.  Another  method 

uses  electroabsorption  (photon-assisted  tunneling)  to  shift 

2 

the  absorption  edge  to  longer  wavelengths  .  This  has  been 

3  4 

employed  in  reverse-biased  GaAs  and  Ga^  A1  As  devices  . 

A  recently  proposed  scheme  is  to  form  an  alloy  between  two 

materials  with  different  crystal  structures.  In  the  region 

where  the  alloy  undergoes  a  phase  change  from  one  structure 

to  another,  the  absorption  edge  is  lowered  due  to  antisite 

disorder.  This  behavior  has  been  observed  in  (GaAs) ^_xGe2x 

5 

alloys  .  In  the  present  work  we  discuss  and  present  experi¬ 
mental  evidence  for  what  appears  to  be  another  method  of 
obtaining  longer  wavelength  response. 

The  scheme  is  essentially  the  inverse  of  a  photon  emis¬ 
sion  process  previously  described**.  It  employs  the  electric 
field  at  the  interface  of  a  hetero junction  to  tunnel  charge 
carriers  across  an  interfacial  energy  gap  smaller  than  the 
energy  gaps  of  the  semiconductors  on  either  side.  Photo¬ 
response  is  obtained  by  the  photoexcitation  of  electrons 
from  the  valence  band  of  one  material  to  virtual  states. 


and  then  electron  tunneling  into  the  conduction  band  of  the 
other  material.  Conversely,  electtons  from  the  valence  band  of 
one  material  can  tunnel  into  virtual  states  and  then  be  photo- 
excited  into  the  conduction  band  of  the  other  material.  The 
latter  process  is  equivalent  to  hole  tunneling,  and  both  are 
illustrated  schematically  in  Figure  4.1. 

4.1  ABRUPT  HETEROJUNCTION  MODEL 

One  of  the  he tero junction  systems  which  exhibits  the 
required  energy  band  line-up  is  In^_xGaxP-GaAs.  An  abrupt 
he tero junction  model  for  this  system  is  shown  in  Figure  4.1. 
Several  of  the  parameters  of  the  model  vary  from  sample-to- 
sample  and  were  determined  in  the  following  ways:  The  energy 
gap  of  GaAs  was  measured  by  photoconductivity  on  high-purity 
epitaxial  samples  to  be  1.39eV.  The  energy  gap  of  In.  Ga  P 

A 
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was  determined  from  its  dependence  on  lattice  constant  , 
where  the  latter  was  obtained  by  double-crystal  X-ray  dif¬ 
fraction  scans.  On  this  basis  the  energy  gap  of  the  lattice- 
matched  alloy  In^  4gGag  was  1.90eV.  The  separations 
between  the  electrochemical  potential  energy,  c,  and  the 
band  edges  were  determined  from  Hall  measurements:  first 
on  the  GaAs  substrates  and  then  on  In.  Ga  P  layers  grown  on 
semi-insulating  GaAs  by  liquid  phase  epitaxy  under  the  same 
conditions  as  those  grown  on  n+  substrates.  The  built-in 
potential  and  the  band-bending  on  either  side  were  obtained 
from  capacitance-voltage  measurements.  From  these  measure¬ 
ments  all  the  parameters  of  the  model,  including  A  6V,  A  £c, 
and  A£I#were  determined  experimentally. 


n-GaAs 


For  a  lattice-matched  hetero junction  with  a  nominal  {111} 
interface,  the  band  offsets  were  A£v=0.08eV,  A  £c=0.59eV,  and 
A^sl.SleV.  Although  these  results  are  consistently  obtained 
experimentally,  it  should  be  pointed  out  that  they  are  not 
close  to  the  expected  theoretical  values.  Using  a  dielectric 

Q 

model  ,  for  example,  we  calculate  ionization  energies  for  GaAc 
and  In^  .49Ga0.51P  of  5.68eV  and  5.91eV,  respectively.  These 
ionization  energies  and  the  experimental  energy  gaps  yield 
band  offsets  of  A£v=-0.23eV,  A£c=0.28eV,  and  A  equals 
the  bandgap  of  GaAs.  Thus,  below  bandgap  photoresponse  is 
not  predicted  from  this  dielectric  model.  Although  there  is 
probably  no  reason  to  expect  good  agreement  between  theory 
and  experiment,  we  have  some  preliminary  evidence  that  the 
valence  band  offset  is  in  the  predicted  direction  for 

g 

In^^Ga^-GaAs  hetero  junctions  with  nominally  charge-neutral 
{211}  interfaces.  Thus,  at  least  part  of  the  difference  in 
band  lineups  may  be  due  to  dipoles  associated  with  the  re¬ 
sidual  charge  at  {111}  interfaces. 

4 . 2  PHOTOVOLTAIC  RESPONSE 

Figure  4.2  shows  the  room  temperature  photovoltaic 
response  for  a  {111}  hetero junction  at  0  and  -0.5  volt  applied 
bias.  (These  diodes  were  rather  leaky  and  exhibited  sub¬ 
stantially  increased  noise  at  higher  reverse  voltages) .  The 
salient  features  are  a  below-bandgap  peak  with  a  half  power 
point  at  1.31eV  which  increases  substantially  with  reverse 
bias  as  expected  for  tunneling.  The  1.31eV  half-power  point 
is  equal  to  the  value  obtained  for  the  interfacial 


Figure  4.2  Relative  photovoltaic  response  as  a  function  of  photon  energy  for 
the  he tero junction  of  Figure  4.1  showing  below  bandgap  response  at 
two  bias  voltages. 


energy  gap,  A  c>1  .  This  exact  agreement,  however,  should  not 

be  taken  seriously.  Also,  in  this  hetero junction  the 
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In1_xGaxP  epitaxial  layer  was  doped  to  produce  2.5x10  cm 

18 

holes,  while  the  GaAs  substrate  had  1.4x10  cm  electrons. 
Thus,  most  of  the  potential  was  in  the  In.  Ga  P  where  it 
should  enhance  the  lower-effective-mass  and  higher  probability 
electron  tunneling.  The  GaAs  absorption  edge  is  apparently 
shifted  to  higher  than  bandgap  energy  due  to  the  high  electron 
concentration^-0 . 

To  test  the  concept  that  the  below-bandgap  response  is 

associated  with  tunneling  across  the  interfacial  energy  gap, 

several  other  experiments  were  performed.  When  the  doping  of 

17 

the  p-type  Inj_xGaxp  was  increased  from  2.5x10  to 
18  —3 

6.0x10  cm  on  the  same  n-type  GaAs  substrates,  the  below- 
bandgap  response  was  greatly  reduced.  Increased  doping  on 
the  p-side  is  expected  to  reduce  electron  tunneling.  Since 
larger-effective-mass  hole  tunneling  from  the  n-side  is  much 
less  probable,  the  smaller  below-bandgap  response  is  consistent 
with  a  tunneling  mechanism.  To  eliminate  the  possibility 
that  the  below-bandgap  response  was  due  to  impurity  tran¬ 
sitions  in  the  space-charge  region  on  either  side  of  the 
interface,  photoconductivity  measurements  were  performed  on 
the  epitaxial  layer  and  on  the  substrate.  These  measurements 
indicated  no  below-bandgap  photoresponse.  Except  for  the 
coincidence  of  the  response  energy  with  a£t  ,  we  cannot. 


however,  rule  out  the  possibility  of  impurity  transitions 
at  the  interface  itself. 

4.3  CONCLUSIONS 

In  conclusion,  we  have  observed  below-bandgap  photovoltaic 
response  in  In^_xGaxP-GaAs  p-n  hetero junctions  which  is 
apparently  due  to  photon-assisted  tunneling  across  the  inter¬ 
facial  energy  gap  between  the  two  materials.  Although  the 
observed  response  for  {111}  interfaces  at  1.31eV  is  not  at  a 
particularly  useful  wavelength,  it  should  be  possible  to  apply 
the  same  idea  for  infrared  detectors  at  longer  wavelengths. 

We  are  currently  investigating  the  possibility  of  obtaining 
smaller  interfacial  energy  gaps  in  the  In^_xGaxP-GaAs  system 
by  examining  interfaces  with  different  orientations.  Other 
lattice-matched  hetero junction  systems,  of  course,  should  also 
be  examined.  Our  results  indicate,  however,  that  hetero¬ 
junctions  with  the  required  band  lineups  can  not  be  predicted 
with  great  reliability. 
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BELOW  BANDGAP  EMISSION  FROM 


ZnSnP„-GaAs  HETEROJUNCTIONS 


Part  of  this  program  involves  the  use  of  the  ordered 
chalcopyrite  compound  ZnSnP2  as  a  template  to  obtain  long- 
range  order  in  the  alloy  Ing  5GaQ  jP.  The  compound  ZnSn?2 
lattice  matches  GaAs  and  these  two  materials,  because  of 
their  different  P  and  As  anions,  produce  a  he tero junction 
lineup  with  interesting  potential  applications^". 

Figure  5.1(a)  shows  the  energy  band  diagram  of  a  ZnSnP2-GaAs 
he tero junction  under  forward  bias.  In  this  configuration 
electrons  in  the  conduction  band  accumulation  region  can 
tunnel  into  the  forbidden  gap  and  recombine  with  holes  in 
the  valence  band  accumulation  region  (and  vice  versa) . 

This  should  produce  infrared  emission  below  the  bandgap  of 
either  material.  Figure  5.1(b)  shows  the  actual  emission 
spectra  from  such  a  structure,  indicating  that  such  a  scheme 
is  feasible.  The  larger  peak  at  higher  energy  is  emission 
from  the  bandgap  of  GaAs,  while  the  smaller  peak  at  lower 
energy  is  apparently  emission  from  the  residual  bandgap 
between  GaAs  and  ZnSnP2.  We  are  currently  investigating 
this  process  in  more  detail. 
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Figure  5.1  (a)  Abrupt  hetero junction  model  for  ZnSnP2-GaAs 

under  forward  bias,  (b)  Emission  spectrum  from 
structure  in  (a) . 
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